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Filial  Roport  Siiiiiiiiary 

High-speed  parallel  processing  has  made  possible  many  real-time  applications  of  digital  signal  proc-ss- 
ing.  Our  research  has  been  mainly  concerned  with  the  design  of  algorithms  within  this  context  We  shall 
summarize  the  problems  studied  and  the  contributions  made  below.  I'ietailed  descriptions  of  our  findings 
will  be  found  in  the  references  listed  in  Part  II  of  this  report. 

In  references  [2],  [3],  and  [6],  adaptive  array  processing  problems  are  studied  with  the  goal  of  finding 
iiumerically  efficient  updates  of  the  singular  value  decomposition  of  the  data  matrix.  This  problem  is  im¬ 
portant  in  the  tracking  of  time-varying  subspaces,  a  generalization  of  the  notion  of  separating  a  time  series 
into  signal  and  noise  components.  The  main  contributions  include  a  numerical  stabilization  technique  which 
eliminates  the  error  buildup  problem  in  a  computationally  efficient  manner  within  the  context  of  practical, 
real-time  computation. 

In  references  [1],  [7].  and  [8],  the  classical  problem  of  spectral  factorization  is  addressed  using  a  .Newton- 
Haphson  algorithm.  By  enforcing  symmetry,  a  special  Euclidian  algorithm  with  .some  remarkabh'  (uoperties 
emerges.  It  involves  updates  which  ar"  formally  equivalent  to  a  backwards  Levinson  algorithm,  which  may 
have  inherent  parallelism.  Furthermore,  the  set  of  minimum  phase  polynomials  is  preserved  under  the 
update. 

In  references  [4],  [o].  [10],  and  [13],  another  classical  problem  is  studied.  A  quite  general  theory  of  what 
is  possible  in  the  bias/variance  tradeoff  in  using  quadratic  e>-:imators  for  the  power  spectrum  is  developed. 
A  very  interesting  time-frequeticy  localization  bound  can  be  found  in  [4].  Iti  [13],  a  new  multiple-window 
quadratic  estimator  with  some  retnarkable  algebraic  properties  is  described.  The  windows  involve  a  form  of 
time  division  multiple.xing  with  circular  shift  constraints.  These  constraints  allow  for  an  algebraic  decom¬ 
position  of  tlie  problem  of  fitiding  orthogonal,  time-shifted  windows  with  maximum  frecpiency  selectivity. 
While  the  estimator  rerjuires  the  computation  and  storage  of  only  one  window  functioti.  its  performanci'  i.-; 
very  simlar  to  that  of  the  Thomson  estimator,  which  uses  the  .^lepian  window  set. 

In  references  [14]  and  [1')].  a  problem  of  designing  wavelets  with  maximum  frequency  sj^lcft  ivity  '^ub  joc) 
to  a  constraint  involving  the  length  of  the  underlying  FIR  filter  and  the  regularity  (or  nutiiber  of  zeros  at 
c  =  -1)  IS  characterized  and  solved.  The  solution  involves  characterizing  extremal  wavelets  in  tiTiiis  of  their 
zeros  on  the  unit  circle  and  tlu’ii  finding  the  zeros  of  tin'  extremal  filter  via  a  generalization  of  the  classical 
problem  of  numerical  quadrature  A  doubly  symmetric  Euclidian  algorithm  for  the  robust  determination  of 
the  magnitude  of  the  wavelet  frequency  function  is  al.so  develojied 

In  referetices  [9],  [11],  atid  [12],  some  parametric  modeling  jiroblems  are  addressed.  In  [9],  a  unified 
approach  to  low-order  modelling  which  involves  projections  within  high-order  tnodels  is  found  to  include 
several  of  the  classical  model  reduction  methods. 
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